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Abstract 

One of the simplest pseudo-Hermitian models with real spectrum (viz., square- well 
on a real interval X of coordinates) is re-examined. A PT— symmetric complex 
deformation C of X is introduced and shown tractable via an innovated approach 
to matching conditions. The result is surprising: an arbitrarily small deformation 
1 ^ C implies a sudden collapse (i.e., the spontaneous PT— symmetry breaking) of 
virtually all the spectrum (i.e., up to its low-energy part). 
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1 Non-Hermit ian Hamiltonians and their spectra 



Thirty five years ago, Bender and Wu [1] published an extremely exciting discov- 
ery that certain bound-state problems may be much better understood when one 
drops the "obhgatory" Hermiticity assumption H — and admits that a coupling 
constant > in Schrodinger equation 

Hig) |V'„(g)) = EM IMg)), n = 0, 1, . . . (1) 

is analytically continued to a complex value g & (D. In this perspective Bender and 
Wu worked, for definiteness, with the quartic anharmonic-oscillator Hamiltonians 

H{g) = H^'\g)=P' + px^ + gx\ (2) 

and demonstrated that the separate (though, in general, complex) spectra {En{g)} 
may all be interpreted as the sets of an intersection of all the Riemann sheets of a 
single analytic function IE'^'^\g) with a corresponding "hne" of a constant g. Many 
years later, similar observations were made and verified for the cubic model 

H{g)^H^'\g)^f + fx' + igx' (3) 

(see the review [2] for more details) etc. In all of these models, the costs of the gen- 
eralization H ^ proved much lower than expected. For all the nonzero couplings 
g 0, all of their complex "exceptional points" (EP) [3] proved well separated in 
the complex plane of g for both m — A and m — 3. 

A new important development of the subject emerged cca six years ago when 
Bender and Boettcher published their letter [4]. Having extended their attention to 
the whole class of the power-law models 

H{g)^H('+'\g)^f + fx' + gx'{ixY, S > 0, g>0 (4) 

[reproducing the above special cases (2) and (3) at 6 = 2 and 6 = 1, respectively], 
they summarized several existing perturbative and numerical experiments (for illus- 
tration one could cite, e.g., [5, 6, 7]), complemented them by numerous new WKB 
arguments and conjectured that after one introduces a suitable Hilbert space, all 
the Hamiltonians (4) (with, for simplicity, / = 0) will possess the purely real (and 
discrete and, from below, bounded, i.e., "observable- like" ) spectra, in spite of their 
manifestly non-Hermitian origin. 
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For the pioneering conjecture of this type, several rigorous proofs {cf., e.g., the 
Fourier-transformation results [7] for eq. (2) at g < 0, or the manifest reality of 
perturbative energies after re-summation for eq. (3) at the real g [5]} were al- 
ready available and many other had only to come {cf., for illustration, proofs in 
the difficult case of eq. (3) in [8]}. Nevertheless, the core of the message deliv- 
ered by Bender and Boettcher lied in the emphasis attributed to the parity-plus- 
complex-conjugation symmetry (conveniently called PT— symmetry) of their sample 
non-Hermitian Hamiltonians H^'^^^\g) with real spectra. This inspired an exten- 
sive subsequent study of the structure of the relationship between the reality of the 
spectrum and the T'T— symmetry of the underlying non-Hermitian Hamiltonian [9] . 

We intend to contribute to the latter effort by the description of an exactly 
solvable example which exhibits a rather counterintuitive enhanced sensitivity to a 
very small change of its coordinate domain. We shall start from an overall review of 
the state of the art in section 2 where we emphasize the theoretical importance as 
well as some practical weaknesses of the pseudo-Hermitian constraint imposed upon 
the non-Hermiticity of the Hamiltonians. 

In the next section 3 we return to the study of quantitative characteristics of 
the specific, differential-equation models where the current and robust property of 
Hermiticity H — is being replaced by the PT— symmetry which may be fragile 
[10, 11]. We restrict our attention to the most elementary non-Hermitian square- 
well model (of rcf. [12], with real spectrum) and extend its scope slightly by the 
replacement of its usual domain (viz., a finite interval 1) by a broken line (or by any 
other smoothly deformed curve C - cf. Figure 1) in the complex plane of coordinates. 

Our main mathematical results are presented in section 4 where our new method 
of solving the matching conditions is shown applicable to an explicit qualitative 
description of the structure of the bound states in the broken-path regime. We 
explain in detail how our "moving-lattice" method decisively facilitates the global 
analysis of the matching conditions. 

Our key physical message is finally formulated in section 5 emphasizing that our 
innovative geometric interpretation of the matching conditions offers the rigorous 
proof that in our example, an arbitrarily small imaginary shift — > i a; of the match- 
ing point causes a non-perturbative, sudden complexification of all the high-energy 
part of the spectrum. In the currently accepted terminology this means that the 
vast majority of the wavefunctions encounters a spontaneous breakdown of their 
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PT— symmetry. We shall conclude in section 6 that this symmetry is manifestly 
fragile in our particular model. 





2 PT— symmetry 
2.1 Variational picture 

In the early stages of study of PT— symmetric quantum mechanics people tried to 
understand the complexified Schrodinger Hamiltonians of the type (4) as models on 
the real line, with a pair of the real and imaginary potentials of a definite behaviour 
with respect to the parity V. Along these lines one arrives at an introduction of the 
two harmonic-oscillator- type bases { jn*^^-*) } (with definite, fixed parities ^^'') and 
transforms the PT— symmetric differential Schrodinger equations H |^) = E 
with a pre-selected normalization of j'^) = (jn^'^^)'ilj^^ + i \n^~^)'ijj^^^ into the 
variational-like real and partitioned matrix problems containing arrays tf^^^ of the 
real wavefunction components il>^\ 

( \ ( \ 

(5) 

The infinite-dimensional submatrices A — and D — are real and symmetric 
but the spectrum itself need not be real at all [13]. These considerations inspired 
Mostafazadeh who conjectured, in a series of papers [14], that the Bender's and 
Boettcher's PT— symmetric quantum mechanics should be classified as a mere special 
case of the more universal pseudo-Hermitian quantum mechanics, the origins and 
foundations of which might be traced back to Dirac et al [15]. In such an overall 
setting, he proposed to weaken the VT symmetry of the Hamiltonians to their mere 
pseudo-Hermiticity 

H^ = riHr]~^, r] = if (6) 
where we may set, in our particular example (5), 

In the light of the well known [16] huge ambiguity of the assignment of the "metric" 
77 to any given pseudo-Hermitian Hamiltonian H ^ H\ Mostafazadeh also proposed 
that the "natural" choices with the indeterminate parity-like metric operators [like 
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rj-p in eq. (7)] should be all replaced by any (i.e., very often, non-diagonal and 
strongly Hamiltonian-dependent) positive definite alternative rj^ > 0. In parallel to 
Mostafazadeh, similar conclusions have been reached in refs. [17] and [18] defining, 
in the present language, the particular positive definite metrics rj^ —VQ and ?7_|_ = 
CV using the additional symmetry generators Q and C of quasi-parity and charge, 
respectively. 

The latter procedure enables us to call al the similar non-Hermitian Hamilto- 
nians H "quasi- Hermitian" because, in accord with the review [16], the positivity 
of 77+ > suppresses many interpretation difficulties and leaves the quasi-norm 
\\(p\\ — (</7|77+|</?) real and non-degenerate. This makes the corresponding Hamil- 
tonians compatible with their standard quantum-mechanical probabilistic tractabil- 
ity [19]. 

2.2 Square-well illustration 

One of the main sources of inspiration for the selection of potentials in Schrodinger 
equations (say, in the Coulombic form) is the principle of correspondence which al- 
lows us to extend and transfer to quantum mechanics the experimental experience 
gained during centuries in the common, macroscopic world. A counterintuitive char- 
acter of many quantum phenomena allows us to search for some new and unusual 
Schrodinger equations, e.g., by a complexification of their axes of coordinates 1? — > C 
[4]. Alternatively, we may obtain manifestly T'T— symmetric equations 

il){x) ^ E il){x) (8) 

by staying on the real line and "deforming" the shapes of V{x) = +V{—x) and 
W{x) = —W{—x). Samples of both these approaches may be found, e.g., in ref. [10] 
(considering mainly the asymptotically power-law potentials) or [20] (paying atten- 
tion to some exponentially confining forces) or [12] (where an even steeper, infinitely 
deep square well has been complexified) or [21] (where the mathematical proper- 
ties have been discussed for the next-step models with delta-functions mimicking 
the "infinitely thin" square wells). The square- well model represents a reasonable 
phenomenological compromise exhibiting, as a bonus, an important merit of exact 
solvability. We shall pick up it in what follows, interpreting the infinitely deep real 



dx^ 



+ V{x)+iW{x) 
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part of the potential 



V{x) = < 



+00 


+00 



for 



X > 1 
-1 <x <1 
X < -1 



(9) 



as a requirement that all the wavefunctions vanish at re = ±1. We break the Hcrmitic- 
ity of our Hamiltonian by adding the imaginary and PT— symmetric finite interaction 
term with coupling Z > 0, 



W{x) 



+iZ, 
-iZ 



for 



Re x < 

Re a; > 



(10) 



Once we assume that the interval of the coordinates x remains purely real, the 
spectrum of energies En — En{Z) proves discrete and real at all Z < Zcruicai ~ 4.48 
(cf. ref. [22]). It smoothly converges towards the well known square- well energy 
levels En(0) = (n + ifn'^/A in the Hermitian limit Z ^ 0. 

In an extension of the above square-well model we shall now assume that the 
interval of the coordinates x is deformed to complex plane. The corresponding gen- 
eralized, T'T— symmetric (i.e., left-right-symmetric) curve C is sampled in Figure 1. 



3 Exact solvability of the new model 

Once we define our potential (10) in the whole complex plane of a; G(Z7, solutions iIj{x) 
will be analytic in both its half-planes. The only distinctive feature of our present 
generalization (—1, 1) = X — > C lies in the requirement of the matching of the left 
and right branches ip^{x) of our full wavefunction iIj{x) at a point Xq = iuu on the 
imaginary axis. This enables us to postulate the matching rules 

ip_{iu) = ip^{iu) = 1, dxip_{iuj) = dxip+{iuj) = iA (11) 

in terms of an auxiliary real parameter A e (—00, 00). 



3.1 Re-parametrization of the matching conditions 

As long as our potentials V and W are constant almost everywhere, the general 
solution of our differential Schrodinger equation (8) may be put equal to a sum of 
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the hyperbolic sine and cosine. The left and right solutions ipifix) are different, 
having to vanish at the different boundary points x ^ ^1, 

V'-(x) = sinhK*(l + 'ijj+(x) ^ sinh k(1 - x) . (12) 

With n = s — it, the values of the two free real parameters s and t will be determined 
by the differentiation in eq. (8), 

E^f-s"^, Z^2st. (13) 

As long as a change of the sign of k would influence just the (arbitrary) sign of the 
overall normalization coefficients R±, we conveniently restrict our attention to the 
quadrant of s > and t > (note that we fixed the sign of Z > in advance) . The 
insertion of the right and left solutions (12) in the matching conditions (11) gives 
the following complex (and transcendental) algebraic equations, 

L sinh K*{1 + iuj) = R sinh k{1 — icu) = 1, 

K*L cosh«;*(l — —kR coshK(l — iu) — iA. 

Their solution is our main task. In the first step, we can get rid of the redundant 
constants by taking the ratios, 

cotanh«;*(l + iu;) — —k cotanh«;(l — iuj) — iA. (14) 

The former equal sign is trivial while the latter one represents a complex equation 
which defines the real parameter A and inter-relates the two real and positive pa- 
rameters s and t in addition. As long as we have Z = 2st, this should determine all 
their admissible values. 

Changing our notation and putting uj = tarup with G — 7r/2,7r/2), let us now 
introduce two auxiliary linear functions S — S{s, t) and T — T{s, t) defined by the 
elementary two-dimensional rotation 





(15) 



cos (p —sm(p 
sin (f cos ip 

where the angle of the rotation measures also the upward shift of the matching point 
in the complex plane of our complex coordinates. In this notation we may re-write 
our matching constraint (14) in the form 

A . , f S T \ ^ , ^ S T \ 
A smh I = If + zs) cosh i (16) 

\COS (fi cos (fi J \COS(fi cos (fi J 
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which admits a facihtatcd separation of its real and imaginary part. The value of A 
drops out of their ratio which may be further re-arranged to represent our matching 
condition in the most compact real form 

/ 2^ \ / 2T \ 

s smh = —t sin . (17) 

\cos (p J \cos (p J 

In the limit — > the latter equation coincides with the cu — prescription of ref. 
[12]. At the generalized u ^ the replacement of s and t hy S and T via eq. (15) 
converts our new and more complicated matching formula (17) into its final form 

CO + qIt) sinh a , , 

T = (7- ^-T — — , u; = tan(/j 18 

1 ~ ^'(''"j ^ smh a 

where we abbreviated a — 2S/ cos</?, r = 2T/ cos</? and q — q{t) — —1/ sinr. This 
equation is an implicit definition of a certain set of curves r = ©(cr) in the a — t 
plane. In principle, the knowledge of these curves would enable us to find all their 
intersections ((7jk, Tk), A; = 0, 1, . . . with our original constraint t — Z/ (2s). 

3.2 The lattice-moving method of solving eq. (18) 

Our present key idea is that the function q — q{t) is periodic, i.e., it remains constant 
on a discrete lattice C of its argument r. In this spirit we shall split the real axis of 
T into intervals of the length 27r numbered by an integer k. Then we introduce the 
second variable p = ±1 marking the right and the left half of each of these intervals, 
respectively. This guarantees that at a fixed p the sign of the sine function remains 
the same and equal to —p. Finally, due to the symmetry of each of the sine-shaped 
curves we split the half-intervals in the quarter-intervals marked by another index 
g = ±l, 

T = (2A; + l)7r + p| + g|e = n^,,) (p, , C = ^tau) e (0, 1) . (19) 

As a consequence, our parameters ^(r) become represented by the functions which 
are independent of k and q, 

^{P, = —. = + /./o^ • (20) 

sm T^k,q){p,0 cos(7r^/2) 

This means that the parameters g remain constant over all the lattices vC(po,€o) 
points T(k,q){po, Co) where the sign po and the parameter are temporarily fixed. 
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3.2.1 Verification: Straight-path solution re-visited 

At a; = and C — X, the use of the hmiting, simphfied version 

T — q{t) a smh. a , u — Q (21) 

of our matching condition (18) leads to an enormous simplification of the construc- 
tion performed in ref. [12]. There, severe difficulties originated from a strong and 
pronounced r— dependence of the factor q = ^(r) which is a very quickly changing 
function of its argument r. In our present setting, the discretization (19) enables us 
to fix the value of g = Q hj reducing our attention from all the vahics of r to their 
lattices C{^p^^). Treating them separately, one at a time, we only have to keep in mind 
the overaU range of our real constants VL{p,^) = p\^{p,0\ = P^(+1)C) ^ 
This enables us to re-parametrize the matching condition (21), 

= = ^sinhcr, p,C = fixed. (22) 

In the new language, the graph of the function ft a sinh o" is a parabolic curve which 
is oriented up or down at the respective p = -|-1 and p = — 1. As long as we are 
interested in the positive r > 0, we may discard p = — 1 and fix r = T(fc_g)(-|-1, > 
and q{t) = ll(p,0 = ^^(+1,0 > 1- The curves r = Q(-M, ^ ^-sinh cr = ©^((j) then 
shrink in proportion to the growth of ^, proceeding from their broadest ^ = version 
(where \ fl\ — 1) via the narrowing parabolic curves until the degenerate single and 
upwards-oriented half-line in the limit ^ — > 1, i.e., fl — > oo. This is illustrated in 
Figure 2 where the unhmited shrinking of the curves is sampled at ^ = 0,0. 5, 0.9 and 
i = 0.99. 

When we zoom out a stripe of r = T(^k,q){p, at a fixed /c = 30 in Figure 2, wc get 
Figure 3. In the new Figure the variations of r arc determined solely by the changes 
of ^ and q which are sampled by a few horizontal lines. As long as the right-hand- 
side function ©5(0") depends on both a and ^, equation (22) will be satisfied only at 
the points of intersection of each particular marked horizontal line with another 
particular, ^—assigned parabolic curve. In this manner the points of intersection 
(cr^^,T^^) in Figure 3 sample the graphical solution of the matching condition (22). 
As long as we choose a fairly large stripe number k = 30, the parabolas of Figure 2 are 
represented by the almost straight and almost vertical lines in Figure 3. This makes 
the identification of all the intersections particularly easy. We see that the points of 
intersection form the horizontally prolate ovals, each of which being confined within 
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its k = ko and p = +1 stripe, and not exceeding the interior of the "maximal", 
|0| = 1 parabola. It is obvious that the horizontal lines (= lattices of r) as well 
as the more or less vertical parabolas move smoothly with the growth of ^. The 
resulting picture reproduces precisely our old graphical proof [12] of the existence of 
solutions at a; = 0. Our present new discretization method appears to offer a feasible 
extension of this proof and analysis to a; 7^ 0. 



4 Bound states at 7^ in graphical representa- 
tion 

Once we wish to determine the spectrum of the square-well energies En = t^ — s^^ at 
any a; 7^ 0, we have to find all the real values of s = s„ and t = tn which satisfy both 
the constant-coupling constraint (13) andihe matching condition (17). In the first 
step, let us re-express the former elementary hyperbolic-curve correlation Z = 2st in 
the new variables a and r. 



4.1 The Z— dependent hyperbolic-curve constraint 

Rotation (15) implies that under the assumption > we have r > while the 
sign of (7 may be both positive and negative. Alternatively, the choice of a; < 
would imply that we must keep a > while the sign of r is allowed to vary. This 
means that one of the two hyperbolas defined by the rule Z = 2st may be discarded 
immediately. Of course, in our innovated notation we must describe these hyperbolas 
by the slightly less transparent rotated quadratic equation 

9 „ cos2(/3 9 AZ ^ 
+ 2t(t - (7^ — = 0. 23) 

sin 2(/9 sin 2(^ cos"' Lp 

At a; = tan 99 > it is easy to select the correct branch defined by the formula 



r = S((7) ^- (uj--\ a+-\ (uj + -\\'' + AX\ = , . (24) 

^ ' 2\ u)) 2 V V cjy ' sin2v5 cosV 

In parallel, at a; = tan = —a; < we must use the different formula 



a = T{T) = U--u;) t + IJL + -)\'^ + AY\ = ■ , 2 ■ (25) 

/ Va; / z V V ouj smZip cos^(/9 

In the other words, we must treat the up and down shifts i a; of the matching point 
separately, reflecting the fact that we already broke the symmetry between the half- 
planes of coordinates x & Chy having chosen the positive coupling Z > in advance. 



9 



4.2 The second, matching constraint 

The application of the lattice-shifting technique of section 3.2 may be extended to 
both the positive and negative cu. The variable r = T(^k,q){p,0 remains represented 
by the same function of the interval selector k, of the two sign- variables q — ±1 and 
p — ±1 and of the continuous ^ varying in the compact interval (0, 1). 

4.2.1 Moving lattices 

Let us now select a; > and keep the two auxiliary variables p and ^ fixed. This 
restricts the range of our variable r to the lattice C = C{p, ^) where the function 
fl = fl{p,^) = ]9/cos(7r^/2) remains constant. This leads to a decisive simplification 
of our matching condition (18), 

^ = Q(P,$)(^) = ^ 1^ ^^^^-J ' p,i = fixed. (26) 
^ ' l — iluj smh a 

A typical graph of the function 6(p^^)((j) at both p = ±1 and at the minimal ^ = 
and/or = 1 is displayed in Figure 4. With respect to the growth of the lattice- 
characterizing parameter ^ from to 1 it is trivial to see from eq. (26) that 

• at p = -1-1 and cr > 0, the right branch of the well-shaped curve ©(+i^^)((t) > 
is bounded by its perpendicular asymptote at (7oo(0 — arcsinh(l/[a;Q(-|-l, ^)]). 
With the growth of ^ and Q, it inadvertently moves to the left and in the hmit 
of ^ — > 1 it coincides with the vertical half-axis Vright (coo(l) = 0); 

• in parallel, the left branch of the same well moves upwards and coincides with 
its diagonal asymptote Ae/t in the same limit, lim^_+i ©(+i_^)((t) = —ojio] 

• at p = —1 and a < —cTooiO < 0) there exists another hyperbolic well sam- 
pled in Figure 4, with definition r = 0(_i^^)(o") and asymptotes -Die/t and 
(^—dependent) Aright- With the growth of ^ this well moves downwards and 
to the right and coincides with the wedge formed by £^/e/t and Vright at = 1. 

This geometric picture has several consequences. The most important one is that at 
the minimal ^ = the curves of Figure 4 contain the initial points of all the ovals of 
the solutions in the manner indicated by the ^ = point in Figure 3 above. With 
the growth of ^ the similar oval-shaped curves are then being formed at any a;. 
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4.2.2 Four families of half-ovals 

In a continuing description of the structure of solutions of eq. (26) we must distin- 
guish between the positive and negative a. For o" > 0, the analysis is simpler since 
the ovals (or rather half-ovals) as sampled in Figure 3 at a; = can solely exist in 
the stripes with ■p — ^l. With the growth of ^ they open their two g = ±1 branches 
to the left until they attain their maximal width and reach the boundaries of their 
stripes on the vertical axis Kig/it in the limit ^ — > 1. 

At (7 < we have to parallel the above half-ovals by their p — +1 partners which 
start to open to the right at the leftmost curve with ^ = 0. They end their growth 
at ^ = 1 while touching the boundaries of their p = +1 stripes on the diagonal Dieft- 

In contrast to our above uj = exercise in section 3.2.1, the choice of o" < admits 
the existence of another family of the half-ovals within the -Die/t ~ Vright wedge. Of 
course, they can only exist within the stripes where p = —1 and in the domain of the 
sufficiently large r ~ Tq (i.e., at A; > kminimai) where they can originate on the curve 
r = 0(_i o)(o") > Tq. In this domain they form the two subfamilies again, depending 
on whether they originated on the left or right branch of the ^ = curve. 

4.2.3 Two patterns of gluing the half-ovals 

With the growth of ^, the left half-ovals within the wedge -D^e/t — Kia^t open to the 
left, ending their growth at ^ = 1 in the intersections of the boundaries of their 
p — —1 stripes with the left diagonal straight line Di^ft- At these points these half- 
ovals meet their p — +1 partners so that in contrast to the a; = pattern (with a 
series of the separated and closed ovals - cf. their picture in ref. [12]), the resulting 
locus of the solutions forms a wavy, sine-like-shaped hne which oscillates to the left 
and right and moves up to the left along the diagonal Dieft- As long as this curve 
remains confined between its two envelopes ©(±i o)(o'), the asymptotic decrease of 
the amplitude of this wobbling is exponential. A schematic example of such a wavy 
curve appears in Figure 5. 

In the same range of the sufficiently large r, the second, similar wavy pattern 
is formed along the axis Vright- In exactly the same manner it results from the 
gluing of the right p = —1 half-ovals which open to the right and reach the line 
Vright where they find a continuation in the above-mentioned p = +1 half-ovals at 
(7 > 0. In contradistinction to the previous case, the amplitude of the wobbling is 
asymptotically constant. Still, this fact alone is sufficient to exclude this branch from 
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further consideration because the Z = 2st constraint is asymptotically a hyperbola 
with asymptotes at the angles ip = arctana; and ip' = arctana; — 7r/2 with respect to 
the axis Kight- 

In all the remaining domain of the not too large values of r, just smooth per- 
turbations ocur of the a; = pattern of disconnected ovals. At uj > the height 
of the ovals exceeds the height of a single stripe. This is consistent with the fact 
that an inner part of the ovals lies within the Di^ft — Vright wedge and must belong, 
therefore, to a, p — —1 stripe. This also does not contradict to the steady decrease 
of the minimum of the graph of the curve ©(_i^^)((7) since with the growth of ^ the 
solutions of eq. (26) start to exist in the lower and lower p — —1 stripes within the 
wedge. In this way, the resulting loci of solutions of eq. (26) are allowed to form the 
separate ovals in a fully consistent manner, indeed. 



5 Energies 

The ultimate goal of our considerations is achieved. We clarified that an optimal 
strategy of the determination of all the parameters s = s„ and i = i„ in the bound- 
state formula En = — sf^, n = 0, 1, . . . is based on a suitable change of variables 
(s",t) (o", t) which merely re-scales and rotates the original hyperbolic constraint 
Z = 2st and re-expresses all the real deformed-path square-well energies by the 
"rotated" formula 

E = En{Z, ^) = ^ [(r^ - (Jr^ cos 2(p-2 (7„r„ sin 2(p\ cos^ (p . (27) 

This leads to a vital simplification of the matching of wavef unctions. In the real 
a — T plane, the construction of all the physical bound states (if any) is reduced to 
an identification of all the admissible parameters ((7„, t„) with all the intersections of 
a certain pair of curves. One of them is the elementary Z— dependent hyperbola (the 
smooth curve in Figure 5). A sufficiently transparent graphical representation of the 
shape of the second one is more difficult and required in fact the greater portion of 
our previous text. This curve is sampled by its quickly oscillating asymptotic part 
in Figure 5. 
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5.1 A comment on asymmetry between u; > and a; < 

We mainly paid attention to the positive values of the shift a; > pertaining to the 
generic form of the family of the hyperbolae given by eq. (24). They are sampled 
by the smoother curve in Figure 5. The Figure also illustrates a generic pattern 
of the intersection of these hyperbolae with half-oval families confined within areas 
specified by their envelope curves exemplified in Figure 4. 

We did not notice in section 4.2.1 that after reflection of Figure 4 with respect 
to the origin of coordinates a and r, its p ~ —1 and p — +1 envelope curves are 
mapped upon each other. This simplifies marginally the construction and follows 
from the invariance of eq. (26) on the lattices since the simultaneous replacements 
T —T and a —o arc equivalent to Q — > — ^2 while the latter change of sign 
merely means that we have to transform p —p. 

We did not deduce, ibidem, that another simultaneous sign-change of a — > —a and 
u —uj preserves the form of the original, lattice-independent matching condition 
(18). This is more important because at the negative u; = — cD < we would be forced 
to replace the most complicated pattern of Figure 4 (where we always employ the 
positive = |ci;|) by its left-right-refiected copy complemented by the corresponding 
correct hyperbolic branch of curve Z = 2st in its alternative form (25). After the 
left-right mirroring transform it enables us to simplify the situation by returning to 
the original Figure 4 complemented by the trivially modified reflected hyperbola 

a = E{T) = l (i-Co) r-lJ(u + i)\' + 4Y\ ^ _-±l^ . (28) 

^ ^ 2 \u; J 2\\ LoJ ' sin2(^cosV 

Hence, all what we have derived at the positive cu may immediately be transferred to 
the case where cu is negative, without changing the half-oval curves and with the mere 
addition of the second branch (28) of the hyperbola. In Figure 5 this would just mean 
a replacement of the upper hyperbola by its minus-sign partner. Of course, such an 
extension of the whole picture is essentially trivial and it need not be discussed 
separately at all. 

5.2 The breakdown of VT symmetry at high energies 

Our construction of a closed form of the bound states is transparent and, undoubt- 
edly, potentially useful. For the straight path C with cu — (p = and for all the 
values of Z > which are not too large, the square-well model already found in- 
teresting applications in the study of the spontaneous PT— symmetry breaking at 
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the sufficiently large Z [22]. An even more important role of this model seems to 
have emerged within the supersymmetric quantum mechanics [23]. In all these and 
similar applications, our present results simply mean that all the changes caused by 
a shift of a small size [a;[ remain smooth if an only if we do not move to the very 
high energies. 

In contrast to that, an introduction of any non- vanishing shift lo changes the high- 
energy region completely and abruptly. In place of infinitely many real and positive 
energies n = 0, 1, . . . which formed the complete spectrum at a; = 0, the choice 

of any cu — tanc/? ^ makes the number of the real intersections (a"„, r„) finite, 
n — 0, 1, . . . , n^aa;(Z, (/?) with a certain maximal real energy at nmaxiZ,^) < oo. 
The mathematical foundation of this conclusion is almost trivial: Up to a finite 
number of exceptions, the energies may only be generated by the intersections in the 
domain of the large r ^ 1 where both the Z— dependent hyperbolas with a; > and 
uj — —\uj\ <Q have almost the same asymptotic representation. 



1 

a; 











+ 1/ 


\uj\)a 



- + 0((7-^), signa; = ±l, cr < -1 . (29) 



This means that both of them share the dominant term (representing just the straight 
line of their common asymptote Dig, ft) and approach this asymptote at an inverse- 
power rate from above or below, respectively (the former case is illustrated in Figure 5 
displaying just the deviation from the asymptote). 

The same asymptote D^e/t is further shared by both the upper and lower envelopes 
©(=Fi,o)(c) of the second, quickly wobbling curve. Nevertheless, from definition (26) 
we easily derive their leading-order asymptotic form 

T = a ^ +y 1^1 ^ Q (sinh-^ a) , sign u; = ±1 , cr < -1 . (30) 
[a; I smhcr ^ ^ 

This implies that the quick, exponential decrease of both the envelopes in eq. (30) 
guarantees that the wobbling line cannot have any real intersections with neither 
of the two hyperbolic Z— dependent curves (29) with their too slow, power-law rate 
of approach to the asymptote. This is illustrated in Figure 5 as a key message of 
the whole construction and implies that the number of the real energies remains 
finite at any non-vanishing 7^ and Z. In the other words, infinitely many real 
energies which existed at 0; = become "lost" and "dissolved" in complex conjugate 
pairs. This occurs precisely at the moment when (say, in Figure 5) the intersecting 
Z— dependent hyperbola moves (say, due to a slight increase of Z) to the top of 
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a particular half-oval (there, the two energies merge at a "Bender- Wu singularity" 
[1] or "exceptional point" [3]) and, in the next stage, separates from the half-oval 
completely (ref. [22] studied this type of a pairwise complexification of the square- 
well energies at a; = in more detail). 

The disappearance of the real intersections of the two curves in Figure 5 occurs at 
any Z 7^ and implies that the VT symmetry of our wavefunctions becomes broken 
at all the sufficiently large energies. In the other words, our initial choice of the 
form of the wavefunctions does not suddenly represent all the possible bound-state 
solutions. In a way discussed in full detail in our previous a; = study [22] this means 
that all the "missing" bound states must be sought in a certain more-parametric and 
manifestly PT— symmetry-breaking form. 

The only exception in encountered a,t Z — where the message offered by Figure 5 
is different because in the limit Z ^ the upper hyperbolic curve moves down and 
coincides strictly with the horizontal axis. This forces us to return to the very origin 
of our present construction and repeat all its steps under the new explicit postulate 
that t — 0. In this case, the kind reader may easily verify that the Z — result 
proves in fact independent of the value of cu so that all the repeated Z — and cu ^ 
(i.e., non-Hermitian though still T'T— symmetric) construction returns us back to 
the energies which coincide with the well known Hermitian square-well spectrum. 

6 Summary 

All our results are summarized in Table 1 which may be read, first of all, as an 
advertisement of our almost involuntary discovery of an extremely elementary and 
transparent new PT— symmetric model with real energies (cf. the last line). On 
a more general level, the main item in the review Table 1 (viz., its last but one 
line) warns against all the non-critical intuition which might prove misleading in 
the realm of VT symmetric models. In this sense, our results may be perceived as a 
non-numerical complement to numerical experiments of paper [10] where several "not 
entirely smooth" potentials clearly inclined towards a spontaneous T'T— symmetry 
breakdown at high energies. 

On this background we believe that in the nearest future, attention will be re- 
attracted to the real role of non-analyticity in the VT symmetric potentials and 
models, with inspiration by our present rigorous proof that any non-vanishing shift 
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of at Z 7^ makes the PT— symmetry of our square- well model suddenly to break 
down. This breakdown involves infinitely many levels at once, i.e., it occurs in a 
way which seems characteristic for virtually all the exactly solvable analyfic models 
[11, 24]. At the same time, the discontinuity of the breakdown might reflect its non- 
analytic origin, contrasting with the robust survival of the reality of spectra under 
path-deformations in many not too strongly singular analytic potentials [25]. 

We have seen that the square-well model is exceptional in representing a solvable 
laboratory which seems to lie on a very boundary between "robust" and "fragile" 
models with T'T— symmetry. In this sense, our present key message is encouraging 
since the geometric language of our innovated "moving-lattice" method proved ex- 
tremely efficient and seems productive. Its former aspect becomes clear when we 
compare the a; = discussion here and in ref . [12] , while its second property is still 
to be verified in the future. 
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Figure captions 

Figure 1. Integration path C with an upweirds deformation a; > 

Figure 2. The lattice-dependent curve (22) with p = -|-1 at a few ^ 

Figure 3. The oval-shaped matching constraint r = ©(cr) in its ^—discreti- 
zation at a; = and A; = 30 

Figure 4. Two lattice-dependent functions r = 6(p o)(<7) at a; = 0.06 
Figure 5. Damped oscillations of ©(cr) vs. constraint 2st — Z 

Table captions 

Table 1. Cardinalities A^(^^"') and iv(co'^p'ex) ^^ve squcire-well energies E 
with Im£^ = and Im£^ 7^ 0, respectively 
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Table 1: Cardinalities A^(^^"') and ]\[(.'^<^p'-<^^) of the square- well energies E with 
Im £■ = and Im £■ 0, respectively 
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j\j'{complex) 
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OO 
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PT— symmetry broken [22] 
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Figure 1 . Integration path C 
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Figure 2. The lattice-dependent curve (24) 
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Figure 3. The oval-shaped matching constraint 
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